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Abstract 

It is commonly believed that most of the stars born in associations decaying with 
characteristic velocities of stars ~ 10 km/s. For dwarf galaxies the decay can lead to 
ejection of stars from the galaxy. The effect is studied for spheroidal and disk dwarf 
galaxies, and is shown to have substantional observational consequences for disk galaxies 
with escape velocities up to 20 km/s, or dynamical masses up to 10^ Mq. The ejection 
of stars can (i) reduce the abundances of the products of Type la supernovae and, to a 
lesser degree. Type II supernovae, in disk stars, (ii) chemically enrich the galactic halo 
and intergalactic medium, (iii) lead to the loss of 50% of the stellar mass in galaxies with 
masses ~ IQ^ Mq and the loss of all stars in systems with masses IO^Mq, (iv) increase 
the mass-to-luminosity ratio of the galaxies. 

Subject headings: ISM: abundances, galaxies: abundances, galaxies: dwarf, galaxies: 
evolution, intergalactic medium 

1. INTRODUCTION 

Mass exchange between a galaxy and the intergalactic medium(IGM) can influence both the 
chemical composition of the galactic gas a nd the morphology of the galaxy. Several mechanisms for 



gas loss by a galaxy can be distinguished (IShustov et al.lll997l ): galactic wind induced by numerous 



supernova explosions, ram pressure exerted by the IGM, the tidal influence of other galaxies of the 
group, evaporation of gas due to interactions with the hot IGM, and blowing of the gas out of the 
galaxy by stellar radiation. Possible ways for the loss of stellar mass include tidal interactions, the 
ejection of stars due to the statistical mechanism, and the decay of stellar associations. Let us 
consider some of these mechanisms in more detail. 

There have been many studies of the effects of nume rous supernova explo s ions, such as galactic 



fountains, super-bubbles, and winds, on disk galaxies (see lShustov et al.l (119971 ) : ICooper et al.l (120081 ) 



and references therein). The efficiency of these processes for gas ejection depends strongly on 
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the distribution of the gas. Thus, models of galaxies with a stratified interstellar medium (ISM) 
displa y a much higher ni ass-loss efficiency than do models with a continuous distribution of their 
ISM (jCooper et al.ll2008l ). Computations of models with a continuous gas distri bution, in turn, 
provide different results de pending on the distributi o n law (|Mac Low et al.lll989l ). It was shown 



in the theoretical study of iDe Young and Gallagherl (jl99d ) that the fraction of expelled gas in a 
1.4 X 10^ Mq galaxy is ~ 0. 6, but, as the author s note , the presence of dark matter was not taken 
into account. According to llgumenshchev et al.l ( 19901 ) . galaxies with masses exceeding 10^^ M( 



o 



do not have w i nds, a nd, hence, do not lose gas via this mechanism. Based on these computations, 
Shustov et al.l (jl997l ) used a simple approximation for the relation between the mass of the galaxy 
and the fraction of expelled matter in their models: 

Mg 



/esc = 2.4- 0.2 Ig 



M, 



1 







In this approximation, the efficiency of gas ejection becomes unity for galaxies with masses of 
10^ Mq, irrespective of their morphology. According to this model, galaxies with such masses 
should not contain gas. At the same time, gas is almost completely absent only in spheroidal 
and elliptical dwarf gal axies, whereas it can constitute a substantial fraction of the masses of dis k 



and irregular galaxies (IBegum and Chengalur 



2004 



Karachentsev et al. 



2004 



Begum et al.l 120081 ). 



On the other hand, observations have not revealed gas ou tflows into the IGM from galaxies with 



dynamical masses of ~ 10^ Mq (Ivan Evmeren et al.ll2009l ). Thus, the question of the efficiency of 



gas ejection due to galactic winds remains open. 

Tidal interaction may be responsible not only for mass exchange between galaxies during colli- 
sions or close fly-bys and betwe en galaxi e s and the IGM, but also for changes in galactic morphology. 
According to the estimates of iTutukovl (120061 ). every galaxy in a cluster experiences a collision at 
least once during its lifetime. In these collisions, the galaxies may merge, lose their gaseous compo- 
nents, or be disrupted completely. A new galaxy may also form from gas lost in galaxy collisions. 

The ram pressure of the IGM gas, evaporation of gas, and sweeping-out of dust are less effi- 
cient galactic mass-loss mechanisms, though they influence the chemical evolution of galaxies and 
enrichment of the IGM. 

The essence of the statistical mechanism is that, in the case of an equilibrium distribution 
of the stars in the gravitational potential of the galaxy, there will always be stars with velocities 
exceeding the escape velocity. As these stars leave the potential well, the system relaxes to a new 
equilibrium state. Howev er, the timescale for the stat istical mechanism is very large — close to a 
hundred relaxation times (IBinney and Tremainelll987l ). where the latter is 

O.liV , , 

Tdyn , (2) 



Trelax 



InA^ 



where is the number of stars in the system and r^yn the dynamical timescale of the system. Typical 
galactic dynamical timescales are ~ 10^ — 10*^ yr. Even for ~ 10^ Mp) galaxies, the relaxation time 
exceeds the Hubble time. Other mass-loss mechanisms rewiewed by iBinnev and Tremaind () 19871 ) 
for collisionless stellar systems are even less efficient. 
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It is commonly beli e ved t hat most stars are born in associations (see however, the paper by 
Elmegreen and Efremovl (Il996l )). The lifetimes of OB associations from birth to decay is short, of 
the order of several million years. Typical veloc ities of stars acquired during th e decay are of the 
order of 10 km/s, according to various studies (jGvaramadze and BomansI 120081) and obser vations 
(jGiesI 119871 ) . Other estimates limit the velocity range to 2 — 8 km/s (IBrqwn et al.l 119971 ) . The 
virial velocities in low-mass gal axies can be several km/ s (IKarachentsev et al.ll2004r) . and the escape 



velocity lower than 20 km/s (IBovill and Ricotti 



2009 



Dijkstra et al 



200 



I 



In the case of disk 

galaxies, the ordered motions of the galactic matter may facilitate the ejection of stars. The aim of 
our current study is to estimate this effect and observational manifestations in dwarf galaxies. 

In Section 1, we compute the probability of ejection of stars from spheroidal and disk galaxies. 
In Section 2, we present the results of modeling the evolution of dwarf disk galaxies taking into 
account the ejection of stars. Section 3 discusses our results. 



STELLAR EJECTION MECHANISM 



To escape its galaxy, the kinetic energy of a star must be sufficient to bring about the work 
against the gravitational field 

(v + 



(3) 



where v is the instantaneous velocity of the OB association, u the stellar velocity relative to the 
association, and $ the gravitational potential at the location of the association, with the potential 
at infinity being zero. The dynamical timescales for associations 



~ 10^ yr) exceed their lifetimes 



10^ yr) due to the low densities of associations, which are 0.1 Mq/pc^ ( Brown et al.l 119971 ) 



This means that energy equipartition for stars of different masses does not have time to become 
established in the association. Therefore, we will assume that the velocities of stars of any mass 
have isotropic Gaussian distributions with dispersion (Tq^; the probability that a star moving away 
from the center of the OB association will overcome the potential of its galaxy will then be 



x(v,-$) 



(y±H)i>_$ 

2 ^ ^ 



(2vra23)3/2 



exp 



(4) 



The large-scale motion of the ISM in the galaxy increases the fraction of ejected stars due to a 
sort of "slingshot effect". If the ISM participates in the Keplerian motion of the galaxy with circular 
velocity the velocity of a star after decay will be summed with the velocity of the association in 
the galaxy. As a result, the probability of ejection of the star becomes (see the Appendix) 



-5/2 
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(5) 
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where 

V I r , 2$ 

^OB V g Or (j2 3 

Here, the velocity of the association is assumed to be equal to the local velocity of the ISM. The 
reason for this is that the lifetimes of associations are much shorter than the galactic dynamical 
timescale so the velocity of the association does not change appreciably before the association 
decays. The parameter ij in Eq. |5] describes the relative velocity of the Keplerian motion in units 
of o"oB- Expression [5] is also valid for galaxies in which large-scale motions of gas are absent, such 
as spheroidal galaxies, but we must then adopt rj = 0. 

The value of x depends on the mass distribution in the galaxy and the motion of the galactic 
gas. It is interesting to estimate the fraction of ejected stars for some typical configurations of 
galaxies. We will calculate this estimate as an average over the volume of the galaxy weighted by 
their local star-formation rate 

Let us give estimates for a spheroidal galaxy with a Plummer density profile and a disk galaxy with 
an exponential density profile. 

The distributions of the density and potential in an isotropic Plummer sphere depend on two 
parameters — the mass M and the characteristic scale a: 

p = 1^ (1 + rVa')~^/' , $ = (1 + r'/a'yy^ . (8) 



Let us take the model of lFirmani and Tutukovl (jl992h for the volume star-formation rate: 



^ oc . (9) 

The distribution of the dimensionless potential ip over the scale r/a depends only on the parameter 
P, which is defined as the ratio of the typical virial velocity in the galaxy and the velocity dispersion 
in the decaying association: 



/ 2/32 GM 



(l + r2/a2)i/2 ' 



As we noted above, we assume rj = for a spheroidal galaxy. Inserting these relations into Eq. [7] 
and integrating over the volume, we obtain the coefficient x as a function of f3 (Fig. H] left panel, 
solid curve). 

The distribution of the surface density of an exponential disk also depends on the total mass 
and spatial scale: 

" ^ e-^/<^ . (11) 



27ra2 
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Fig. 1. — Left panel: average probability of ejection of stars from galaxies as a function of (3, for disk galaxies 
with an exponential density profile and circular Keplerian rotation (solid) and non-rotating spheroidal galaxies with 
a Plummer density profile (dashed). The central panel shows /3 as a function of the mass o f the galaxy (up to 
10^0). The mass distribution for the galaxies was taken from the catalog of nearby galaxies of lKarachentsev et al.l 
(200J). M25 is the dynamical mass within the 25™ isophote. The parameter j3 was computed as the ratio of the 
rotational velocity of neutral hydrogen (denoted in the catalog Vm) and the velocity dispersion in the decaying OB 
association ctob- Two limiting cases were considered for the latter: ctob = 2 km/s (triangles) and ctob = 8 km/s 
(circles). The solid lines show approximations for /?(M25) for these limiting cases. The dashed lines show analogous 
approximations for the ejection of stars outside the dark- matter halo (see the text). Right panel: dependence of the 
average probability of ejection of stars on the galaxy mass M25. The lines correspond to the same cases as in central 
panel. 



Let US apply an approximate expression for the gravitational potential, with the mass distribution 
taken to be spherically symmetric: 



GM 1 - e"^/*^ 
a r/a 



(12) 



In this approximation, as in the case of a Plummer sphere, the distributions for r/ and ijj depend 
only on /3, which has the same definition: 



7? = /3. 



1 - (1 +r/a) e-"/" 



r/a 



1 - e-""/" 
r/a 



(13) 



Finally, let us take the surface density of the starformation rate according to the SchmidtBl^YKen- 
nicutt model (lKennicuttlll998l l: 

^! oc . (14) 



The function x(/3) for the two cases is plotted in Fig. [T] (left). For disk galaxies, the values 
of /3 for which the ejection of stars is possible are bounded from above by 4 — 5, while this limit 
i s 2 — 2.5 for spheroidal g alaxies. We used data from the Karachentsev catalog of nearby galaxies 
(IKarachentsev et al.ll2004[ ) for galaxies with masses below 10^ Mq to derive the dependence of j3 on 
the galactic mass. For each galaxy with a given mass, the value of /3 was calculated as the ratio of 
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the rotational velocity of neutral hydrogen (denoted in the catalog Vm) and the velocity dispersion of 
stars in the decaying OB association cob- Since cob can be in the range 2 — 8 km/s, the parameter 
/3 is also not determined precisely. The approximate mass dependence of this parameter is 



/3 = 0.025 



M' 



25 



Ma 







0.36 



/ O'OB 
\km/s 



(15) 



In this approximation, we use the dynamical masses of the galaxies inside the 25 mag arcsec"^ 
isophote, denoted M25. Using only M25 for the mass underestimates the contribution of dark 
matter in the outer regions of the galaxy but the approximation [15] remains valid even if we take 
into account the dark haloes, viz. if we increase M25 by factor equals to 1 + i^dm/^h ~ 6, then the 
coefficient on the right-hand side becomes 0.048. The function [T5] for the limiting values of aoB is 
shown in Fig. [1] (middle panel) for the cases when the dark-matter halo is and is not taken into 
account. 

Below, we present results of our numerical modeling of dwarf galaxies in a single-zone ap- 
proximation, taking into account the ejection of stars. We traced the chemical evolution of the 
galactic gas using the abun dances of iron, oxy gen and heavy eleme n ts. We used the n iodel of 
Firmani and Tutukovl (jl992h and the studies of IShustov et al.l (jl997l ) ; IWiebe et al.l (jl998l ) as the 
basis for the one-zone model. 



3. NUMERICAL MODELING 



3.1. Single-zone Galactic Model 



In a decaying association, only stars with lifetimes < 10'' yr end their lives in the galaxy as Type 
II su per novae (SN II). These stars have masses of about 13 Mq and higher (iTutukov and Krugel 
198d ). Some fraction of lower-mass stars, whose number is determined by x, leaves the galaxy 
and is not able to participate in the enrichment of the ISM. Ho wever, this effect may be small, 
since stars with masses below 8 Mq do not explode as SN II (iHashimoto et al.l Il993[ ). so that 
their disappearance influences the chemical composition of the ISM only more weakly. On the 
other hand. Type la supernoy ae (SN la) have delay times until their explosions of 10* — 10^ yr 
(ITutukov and YungelsonI l2002l ) . and their ejecta can strongly influence the enrichment of the ISM 
in metals. The mass fraction of all stars born as a single stellar population with a given initial mass 
function (IMF) <I>(m) and leaving the galaxy is y fl^^^'^ dm^; for Salpeter IMF with index —2.35 
and stellar mass range 0.1 Mq — 100 Mq this fraction is 0.88x- 

In the single-zone approximation, the total rate of stellar mass loss due to the ejection of stars 
for a galaxy with a time-dependent star formation rate ^(t) will be 



M:^{t) = x^it-T,{13MQ)) 



ISA/q 



dm $ 



(16) 
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where Tg = Ts{m) is the lifetime of a star with mass m. 



The ejection of stars also affects the amount of mass participating in star formation and chemical 
enrichment of the ISM. The rate of gas return is 



M^{t) = I dm^it - Ts) (l - ^ ) [1 - xO{IUIq - m)] , (17) 



where = m^{m) is the mass of the stellar remnant, 9{x) = for x ^ 0, and 6{x) = 1 for x > 0. 
The rate of enrichment of the ISM in an element X is determined by the contributions from SN I 
and SN II: 



10" 







+ 



*(i-TSNIa)Pr'^ [1-X] + 

;n.ax r- / 'TYI \ 

dmm{t-T,) (1-— )X(i-r,)+P|^"(m,Z(t-r,)) x 

m 

x[l-xe{13MQ-m)], (18) 

where rsNia is the average delay time between the formation of a binary and the SN la explosion, 
pSNia jg |.]^g mass of element X ejected per explosion, and P^^^{m, Z) is the mass of element X 
produced by SN II with pre-supernova mass m and initial heavyelement abundance Z. 

In this model, we assu med that SN la are produced b y mergers of degenerate dwarfs in close 
binaries. As was shown by iTutukov and YungelsonI (119941 ) . the mean delay time between the for- 



mation of a binary and the SN la explosion is rsNla ~ 10^ yr. The rate of SN la explosions per 
unit star-formation rate is obtained by norr nalizing to the modern SN la rate of 0.003 per yr and 
modern star formation rate of 3 Mq per yr (ITutukov and Yungelsonlll994l ). We assumed that each 
SN la produces P§^^^ = 0.6 Mq of iron (jTsuiimoto et al.lll995l ): we did not consider the production 
of other elements by these supernovae. The yields of SN II as functi ons of the mass of the star and 
the initial heavy-element abundance were taken from iMaederl (jl992l ). 



The balance of the total mass of galaxy, the mass of gas, and the masses of various elements is 
determined not only by star formation, but also by interaction of the galaxy with the ISM; i.e., via 
the galactic wind, dust ejection, and the accretion of intergalactic gas. The se factors were studied 
previously for the same single-zone model (see, for instance, the study of IShustov et al.l (jl997l )). 
However, accretion, dust ejection, and galactic wind were not taken into account in the present 
study [the last factor due to obvious shortcomings of the simple model for the galactic wind, Eq. 
[1]. As a result, the mass balances of various components of a galaxy was given by the equations: 



Mtot 



+ M: 



fb 



(19) 



(1997 


); 


Wiebe et al. 


(1998) 
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3.2. Results of the Computations 

We can see from the dependence x(M25) (Fig. [!)) that the decay of OB associations essentiaUy 
does not lead to mass loss from spheroidal galaxies, even those with the lowest masses. For this 
reason, we restricted our numerical analysis to disk galaxies. We computed four series of models for 
the evolution of galaxies with masses from 10^'^ Mq to 10*'^ Mq and various values of o"ob- For 
comparison, we also computed a series of closed models. The radii of the galaxies corresponded to the 
relation M oi E?. To avoid a large influence of the initial conditions on the burst of star formation, 
the semi-thickness of the protogalactic disk was set to 10 kpc in all computations. The remaining 
parameters for the computed series of models are given in the table. In series A, ctq-q = 2 km/s, 
while in series B and C ctq-q = 8 km/s. In series C, we allowed for the influence of the dark-matter 
halo [see the comments concerning Eq. [15]. Thus, the series B computations enabled us to estimate 
the effects of ejecting stars from the disk into the galactic halo, while series C illustrated the effects 
of ejecting stars from the halo into the IGM. 





(Job = 2 km/s 
(series A) 


Cob = 8 km/s 
(series B) 


cqb = 8 km/s 
+dark halo 
(series C) 


IgM/Mo 


R, pc 


/3 
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/3 
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X 


6.5 


79.5 


2.73 


0.21 


0.70 


0.92 


1.38 


0.6 


7 


141.4 


4.14 


0.09 


1.07 


0.743 


2.08 


0.32 


7.5 


251.5 


6.26 


0.012 


1.62 


0.481 


3.16 


0.13 


8 


447 


9.48 


~ 


2.45 


0.21 


4.78 


0.035 


8.5 


795 


14.25 


~ 


3.72 


0.075 


7.23 


~ 



Figure [2] shows the integrated characteristics of the model galaxies at the end of the computa- 
tions. As expected, including the ejection of stars in the series A models did not result in significant 
deviations of the integrated characteristics from those for the closed model. In the series B and 
C models, for the lowest mass galaxies, the ratio of the ejected mass and the dynamical mass can 
range from 1.5 — 2 (for matter ejected from the halo into the IGM) to 6 (for matter ejected from 
the disk into the halo). The amount of ejected mass obtained in the series B models provides some 
idea of the extent to which the morphology of a galactic disk can vary. For a galaxy with dynamical 
mass lO'^ Mq, the stehar mass in the halo exceeds the mass of the disk by a factor of 1.5 — 2; i.e., 
such a galaxy cannot be considered a disk galaxy. Galaxies with final dynamical masses of 10^ Mq 
can have a mass of 4.5 x 10^ Mq in the halo and 1.5 x 10^ Mq in the IGM. The mass-luminosity 
relation shows that, compared to the closed model, the total luminosity of the lowest-mass galaxies 
can be half an order of magnitude lower, and the disk luminosity an order of magnitude lower. The 
mass-to-luminosity ratio itself increases by more than a factor of 2.5 for the low-mass galaxies of 
series B. For series B galaxies with masses below 10^ Mq, the total luminosity of stars ejected into 
the halo is a factor of two to three higher than the disk luminosity (Fig. [21 middle panel in the 
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upper row). This also provides evidence for the transformation of the galactic morphology from 
disk to spheroidal. 
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I 1 1 1 1 1 1 1 1 I 1 1 1 1 1 1 1 1 I \ 1 1 1 1 1 1 1 




5 5.5 6 6.5 7 7.5 8 8.5 9 5 5.5 6 6.5 7 7.5 8 8.5 9 5 5.5 6 6.5 7 7.5 8 8.5 9 




5 5.5 6 6.5 7 7.5 8 8.5 9 5 5.5 6 6.5 7 7.5 8 8.5 9 5 5.5 6 6.5 7 7.5 8 8.5 9 



Fig. 2. — Dependence of integrated parameters on the disk masses (on a logarithmic scale in units of Mq) in 
four series of models. Solid diamonds correspond to the closed galaxy model, open diamonds to the model with 
aoB = 2 km/s (series A), solid triangles to the model with aoB = 8 km/s plus a dark-matter halo (series C), and 
open triangles to the model with ctob = 8 km/s without a dark-matter halo (series B). The upper row of plots shows 
the ratio of the ejectedmass and diskmass, the ratio of the luminosities of the ejected stars and the disk, and the 
mass-to- luminosity ratio for the disk (in solar units). The panels in the middle and lower rows show the abundances 
of various elements in the gas and stars (averaged over the stellar population). 

Since we did not consider galactic wind and accretion in our models, the star-formation history 
in the galaxies consisted of a single burst which depleted most of the gas: in all the galaxies, the 
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mass fraction of the gas at the end of the computations was 5 — 10%. The relative abundances 
of elements in the gas decrease weakly with the mass of the galaxy and with increasing aoB (Fig. 
121 middle row of plots). One exception is the models at the low-mass end of series B, for which 
the relative abundances of element grow with the galaxy mass by 0.1 dex, compared to the closed 
models. As the plots show, in all series of models, the mass dependence of the abundances converges 
to a minimum for some mass in all series of models. The reason for this may be a competition of 
two processes determining the abundances in the gas: the enrichment of gas by supernovae and 
the return of gas poor in elements by explosions of new low-mass stars. The influence of the latter 
process may be diminished if a galaxy loses a considerable fraction of its low-mass stars. 

The elemental abundances in the stars systematically decrease with the mass of the galaxy and 
with increasing cjob (Fig. [21 lower row of plots). The ejection of stars influences the iron abundance 
most. In a galaxy with dynamical mass 10^ Mq, the iron abundance can decrease by 0.5 dex in 
the disk stars (series B) and by more than 0.3 dex in the halo stars (series C). Figure |2] shows 
the distribution of the relative numbers of stars over the abundances of iron, oxygen, and heavy 
elements for all the computed series of models. Series B (second column of plots) differs most from 
the closed model. Along with the shift of the distribution toward lower abundances, stars heavily 
enriched in iron appear, which formed in the first few billion years after the SN la explosions. The 
characteristic appearance of the stellar distributions is due to the fact that the galaxies in this model 
experience only one burst of star formation, which occurred in a low-metallicity gas. 



4. DISCUSSION AND CONCLUSIONS 

We have studied the influence of the loss of stellar mass on the evolution of dwarf spheroidal 
and disk galaxies. The decay of OB associations was considered as a possible mass-loss mechanism, 
with the decay enabling some stars to obtain velocities sufficient to escape their galaxy. The decay 
of associations is essentially of no importance for the evolution of spheroidal galaxies. The effect is 
also small for disk galaxies with <tob = 2 km/s. Since a value of ~ 10 km/s is thought to be typical, 
we focused our analysis on models with o"ob = 8 km/s. The results of our analysis are as follows. 

1. During the lifetime of an OB association (~ 10^ yr), the most massive SN II (> 13 Mq) are 
able to enrich the ISM in the products of their explosions. Lower-mass stars that leave their 
galaxy do not contribute to the enrichment of the disk ISM, but instead serve as a source of 
elements for the halo or IGM. The same is true of SNIa. 

2. Disk galaxies that had at the onset of their star formation masses of 3 x 10^ Mq contain 
half of their mass in disk stars and the other half in the halo. The halo luminosities in such 
galaxies exceed the disk luminosities by a factor 1.52. We can thus infer that galaxies with 
masses < 10^ Mq that were initi ally disk galaxies cha. n ge th eir morphology to spheroidal. 



According to the classificaiton of Ide Vaucouleurs et al.l (119911) , spiral galaxies are assigned 



morphological indices T = 4 (see also ICorwin et al.l (119941 )). In the catalog of nearby galaxies 
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Fig. 3. — Distribution of the relative number of stars over the current abundances of iron, oxygen, and heavy 
elements. In the left column, the solid curves correspond to galaxies with masses of IO'^-^Mq in the series A models, 
and the thin dotted curve to the closed model. In the other columns, the solid curves correspond to galaxies with 
masses of IO'^-^Mq, IO'^-^Mq, and IO^'^Mq in the series B models, the dashed curves to galaxies with the same masses 
in the series C models, and the thin dotted curves to the closed model. 



(iKarachentsev et al.ll2004l ). disk galaxies have absolute magnitudes not exceeding —13™ (for 
morphological indices from to 7; i.e., including lenticular and irregular galaxies that are 
closest to disk galaxies); this magnitude corresponds to a luminosity ~ 10'' Lq and a galactic 
mass ~ 10*^ Mq. Lower -mass galaxies are classified as spheroidal and irregular. This is 
confirmed by the computations for the models we have adopted here. 

In systems with masses < 10^ -^o, a large fraction of the stellar mass leaves not only the disk 
of the galaxy, but also the halo (Fig. [U right panel). Thus, if extremely low- mass galaxies 
can form at all, they can lose almost all their stellar population to the IGM after their first 
burst of star formation. As a result, a dark-matter halo enriched in gas should be left. This 
scenario may be important for the problem of missing satellites of the Galaxy. 



4. The ejection of stars increases the mass-toluminosity ratio. For galaxies with total masses 
(disk + halo) of ~ lO''' -Mq, this ratio increases by a factor of 2 — 2.5 (Fig. |2l upper right 
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panel). 

5. The ejection of stars may result in strong variations of elemental abundances in the gas (Fig. 
Omiddle row of plots): along with the systematic decrease of the abundances by ~ 0.05 dex, 
the elemental abundances in the lowest-mass galaxies can increase by 0.1 — 0.15 dex. The 
abundances in stars systematically fall with the galaxy mass decrease, by 0.2 dex. 
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APPENDIX 



We will obtain an expression for the probability of ejection of a star from an OB association 
moving along a Keplerian orbit in a galaxy with a given gravitational potential ^{r). A general 
expression for this probability is 



x(v,-$) 



(2^aoB)3/2 



exp 



2^0B 



(20) 



The velocity vector for a circular Keplerian orbit is 



dr ' 



(21) 



where is a unit vector in the azimuthal direction. The integration domain is in this case described 
by the equation 

ul + ul + 2vu^ + v^ + 2<l>^0 , (22) 

where u± is the modulus of the velocity orthogonal to e^. Solving this equation leads to the 
condition for the velocity component 



€ (-oo,3f?('U^_)] U [R{u^+),+oo) 



U4>± 



-v± 



-2$ - ui 



(23) 



Since the components of the vector u are distributed isotropically and are independent, we can 
integrate the general expression over u^. Introducing the notation 

n.2 



2 ' 
OB 



7] 



0"OB 



2$ 

2 ' 
^OB 



we finally obtain 



-7? - VV; - g 
V2 



erf 



V2 



(24) 



(25) 



If r/ = 0, the general expression 1201 reduces to 



(26) 
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